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Abstract. We analyze the phase space of a particular unified model of dark matter, dark energy, and 
inflation that we recently studied in [1] whose Lagrangian is of the form ^f(X,<j>) = F(X) — V(<j>). 
We show that this model possesses a large set of initial conditions consistent with a successful 
cosmological model in which an inflationary phase is possible, followed by a matter era to end 
with dark energy domination. In order to understand the success of the model, we study the general 
features that unified dark matter (UDM) models should comply and then we analyze some particular 
models and find their constrictions. 
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INTRODUCTION 

The current standard model of cosmology is based, on the one hand, on the existence 
of dark matter as a clustering agent to yield both local galactic dynamics and the large 
scale structure of the Universe and, on the other hand, on dark energy as a substance 
responsible for the cosmic accelerated expansion, and they both comprise around 96% 
of the matter-energy content at present. The rest 4% is mainly in the form of baryons that 
are important to understand what we actually observe. Astonishingly, this cosmological 
model have received very much support from different cosmological probes, such as 
type la supernovae, CMB anisotropics, measurements of Baryon Acoustic Oscillations, 
galactic and cluster dynamics, among others, for a short review see ref. [2]. 

In spite of the above-mentioned success, the standard model of cosmology relies 
on the existence of the dark components that make possible the desired cosmological 
dynamics. So far, we do not know, by certain, the origin of dark matter and dark energy, 
although well-motived candidates exist for their origin. Given this, one is tempted to 
look for alternatives to dark matter and/or dark energy in such a way that the known, 
correct dynamics of the standard model is recovered and, if possible, models having 
smoking guns to be able to discriminate among them with the observations at hand. 
One of the interesting possibilities that has appeared in recent years is that different 
phenomena such as inflation, dark matter, and dark energy could be due to a single 
scalar field [3, 4, 5, 6, 7, 8, 9]. A practical scheme is to look for toy models that can 
accomplish the desired dynamics, and in that sense it is important to understand the key 
elements of the Lagrangian that play a particular cosmic role. Different works [3, 10] 
have analyzed the features that a unified model should have, and have pointed out the 



difficulties to build such a single description of the different phenomena. Recently, a 
unified k-essence model for a particular Lagrangian has been proposed [11] and later 
generalized to a whole class of models in Ref. [1]. It was shown there that these models 
work finely to achieve a cosmological dynamics that emulates that of the standard model 
of cosmology including inflation. However, an initial conditions analysis is missing and 
some key features of the success have not yet been investigated. The aim of the present 
work is therefore to analyze the phase space dynamics of that model to identify the key 
elements to later generalize our results to other unified models, such as unified dark 
matter (UDM) models [12, 13]. 



DYNAMICAL ANALYSIS OF THE UNIFICATION MODEL 

j£? = F(X)-V(0) 



We consider here a scalar field with Lagrangian J?(X,<j)) — F(X) — V(0), where the 
kinetic term F in this generalized class of models is a function of the canonical kinetic 
term X = —g^ v d^d v ^)/2. The cosmological equations of motion for this field in a flat 
Friedmann-Robertson-Walker universe are 



1 



3M2 



(2XF X -F + V) 



and 



j (2XF X -F + V)+ 6HXF X = , 



(1) 



(2) 



where Mpj = 1/SnG and equation (2) is the continuity equation. 

In Refs. [1] and [11] one particular choice is made for these Lagrangians correspond- 
ing to 
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where a = 1 corresponds to de model studied in [1 1]. In those works it is shown that this 
scalar field has the interesting properties to emulate the dark matter and dark energy, and 
in the very early Universe to drive inflation. To obtain this behaviour several elections 
of the constant parameters have to be made, for example for a = 1 then 
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A missing part of the study of these models is the analysis of the phase space of 
its solutions. This can tell us whether the solutions that have these cosmologically 



interesting features are feasible to obtain or not. In other words, it tell us whether the 
initial conditions of the solutions are generic or have to be fine tuned. In Ref. [14] 
a study on the general features of the phase space for the models with Lagrangian 
jf = F(X) — V(0) is presented. Here however, we will carry out a similar analysis 
adapted to the particular choice (3,4). 

The continuity equation, assuming a = 1 becomes 

A(j) + 3HA<j)+m 2 <l) = 3v / 2Aa o sign(0)//, (6) 

where sign(^) is —1 if <j> is negative and +1 otherwise. Performing a change of variable 
to 

z = VA<j) , (7) 
and using the Friedmann equation (1) to substitute the value of H 



we obtain an evolution equation for z in terms of the variables and z. This equation and 
the definition of z together can be treated as a system of first order autonomous equations 
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With the equation of state of the field / written in terms of these variables as 
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It can be seen that the system doesn't have any critical points which means that both 
variables z and are always time dependent. The usual dynamical system analysis based 
on the fixed points is thus not possible in this case. However the system can be solved 
numerically to obtain its phase space, shown in Fig. I. There we have plotted in dotted 
(red) lines those of constant equation of state, the horizontal lines corresponding to 
aty = — l and diagonal lines to ftty = 0. As can be seen the sector of initial conditions 
with big negative values and positive z values evolves towards a solution with equation 
of state near — l which in the phase space corresponds to the left horizontal branch. This 
in the unification models is interpreted as the initial period of inflation in which the 
equation of state of the solution gets close to — l . To see this one can show that equation 
(9) drives z to small values when the constant parameters are in the intervals giving in 
Eq. (5). For the expected value of the field at the beginning of inflation 0, ~ l5Mpj as 
obtained in Ref. [I] under slow -roll conditions, the value of z will be around lO^Mpj 
corresponding to co ~ —0.994. If the system starts in bigger z, equation (9) will acquire 
negative values driving the system to this small value corresponding to a potential 
dominated phase. For more details on the inflationary realization of the model see [l]. 



This solution later crosses the lines corresponding to a equation of state equal to 
(diagonal lines) that in the unification models correspond to the matter domination 
epoch. The time that the system stays in the regime of ftty ~ has to be long in order to 
represent the Dark Matter. This time will depend on the value of the parameters (5) in 
the Lagrangian, and in references [1] and [11] it is shown that these parameters can be 
adjusted in order to obtain this behaviour from a redshift of order 10 10 up until a recent 
time when the transition to aty < has to occur; e.g. Eqs. (5) provide the parameters for 
model a = 1. Finally, the solution evolves towards a second period of aty close to — 1, 
that in the phase space corresponds to the right horizontal branch. The whole behaviour 
occurs also for solutions beginning with big positive values of and negative values of 
Z, in which solutions go from positive to negative values in and live in the z < part 
of the phase space, as can be seen in Fig. 1 . 




FIGURE 1. Phase space for the model, the continuous (blue) lines correspond to the evolution of the 
system. The dotted (red) lines correspond to lines of constant equation of state, both horizontal lines 
corresponding to co^ = — 1 and the four 45° segments correspond to fi)^ = 0. A typical solution represented 
by the thick line approaches first to aty ~ — 1 (inflation), then passes through co^ ~ emulating dark matter 
and finally comes back to aty ~ — 1 at late times as dark energy. 

The analysis of this phase space is important in understanding the dynamics of the 
cosmological solutions of the system (3, 4). We can conclude that an important sector 
of the possible initial conditions can give rise to the behavior needed to unify the 
phenomena of dark matter, dark energy, and inflation. If fact, half of the possible initial 
conditions give rise to the behaviour needed for unification. As was stated in the previous 
paragraph, all the solutions that at early times begin with big negative values of and 
positive values of z, or those that begin with big positive values of and negative values 
of z, achieve a successful unified behaviour. A problem that can be seen in Fig. 1 is the 
crossing of the o)^ = — 1 line, that has been argued in [15] it presents stability problems 
for k-essence scalar fields, however this will occur in a future epoch for the case of 
our Universe, as the current equation of state for the field is expected to be close to 



w = -0.75 [11]. 



Purely kinetic Lagrangian 

The kinetic term in the previous Lagrangian (3) was proposed originally in [16] as 
a purely kinetic Lagrangian, corresponding in our case to the scenario in which the 
potential term (4) is small compared to the kinetic term (3). In this case the system poses 
a shift symmetry — > + 0o which implies that it has only one degree of freedom. This 
simplification makes it possible to obtain an analytical solution to the system [17]. The 
energy density and pressure take the form 

p = (AX) a -M, (12) 
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And defining y = (AX) a ~ 1//2 , this equation can be transformed into: 

% = X«-y), as) 

This differential equation has the analytical solution y = Oo + ce~ 3N , where c is a 
constant of integration. After a few e-folds of expansion the constant term dominates. 
The equation of state can be computed as 

(AX) a - 2aocoVAX + (2a - \)M 
<0 * = (2«-l)pX)«-M] ■ (16) 

From the Eq. (14) we can obtain the critical values of the system as AX\ = and 

AX2 = (X,Q^ 2a , which correspond to systems with equation of state — 1 in both cases. 
To study the stability we expand equation (14) around each critical point. For the first 
one, we expand AX = + £ with £ 1 , then the evolution equation can be approximated 
as 

d(AX) n 6e 
dN ~ l-2a' 

for 1 — 2a > 0, which corresponds to a unstable point. If instead 2a — 1 > then 
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and the critical point is stable for Oo < and unstable for Oq > 0. For example, for 
the case studied in the equations (5) where a = 1 and Oq > 0, the critical point X = is 



unstable. For the cases in which X = is unstable, the value of the late time cosmological 
constant is not p (X = 0) = M, that corresponds to the constant added in the Lagrangian. 

For the other critical point, we expand as AX = a 2 ^ 2a ^ + £, then 

d(AX) 

— — - pa -3e (19) 
dN ' U ] 

corresponding to a stable point. This result is important because the dynamical evolution 
of the system will drive the field to a behaviour similar to a cosmological constant 
at late times, with ftty = — 1 in which the system tends to a value of the field with 

AX2 = CCq^ 2(X ^ corresponding to a density p(X = X2) = C^^ a ^ — M that is a 
combination of constants that yields the dark energy of the model, see [1]. 



CONDITIONS FOR UDM 

In the first section we studied the model (3, 4), which in the late time Universe can 
give rise to the phenomena of dark matter and dark energy. However this model is very 
specific, and therefore the aim of this section is to study the conditions for a more general 
class of scalar fields to reproduce the same dynamical features. If the Lagrangian has a 
general form in terms of the field and the kinetic term Jzf = Jzf (X, 0), its equation of 
state turns out to be 

and the sound speed 
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A sufficient condition for the field to behave as dark matter is that both quantities be 
close to zero [13], leaving to the conditions 
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and 
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There are several Lagrangians that accomplish the above conditions and they have 
been proposed as models for unified dark matter (UDM) meaning that they can behave 
as dark matter and, adding a constant to the Lagrangian, as a combination of dark matter 
and dark energy. To proceed testing different Lagrangians we first consider condition 
(22), to later analyze models with (23). 

An example from the literature is the purely kinetic Lagrangian proposed by Scherrer 
on Ref. [12] corresponding to 



J? = F(X)=F Q +F m (X-X ) 2 . 



(24) 



When the kinetic term is near the minimum (X —Xq)/Xq <C 1, this Lagrangian is 
known to behave as UDM. Another similar example of this type is the ghost condensate 
model [18] that served as an attempt to stabilize k-essence scalar fields when having an 
equation of state (0^ smaller than — 1 . These examples show that the behaviour around 
the minimum is important. Let us analyze the conditions for a general Lagrangian having 
a minimum. In this case it can be expanded as 

J/f (X, 0) = J% + \^i8 2 + ^8 3 + • ■ ■ , (25) 

where S£\ is the zth derivative with respect to X evaluated at the minimum Xq, and 8 is 
the deviation from the minimum, 8 = X — Xq. The constant term Jzfo can be dropped 
from the analysis as can be considered as a cosmological constant. This leave us with 
the condition (22) written as 
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that imposes the condition on 8 /Xq to be small, in other words the deviation from the 
minimum has to be small. The higher order coefficients in the expansion are close to zero 
as long as the first one is, except for very particular cases. In ref. [10] it is concluded 
that for Scherrer's model this deviation 8/Xq has to be smaller than 10~ 16 at the present 
epoch to avoid discrepancies in the structure formation and CMB power spectrum in 
comparison with the observations. 

In ref. [1] it is shown that for the model (3) the deviation 8/Xq is of order 10~ 13 
during the equality epoch and at the present epoch of order 10~ 23 , resulting in a correct 
description of the cosmology. Considering our model in more detail, the pressure and 
density of the purely kinetic part, as given in (12) and (13), turn out to be around the 
minimum AXq = a 2 /( 2a_1 ) as 

s 4 2 ryry (2«-4)/(2a-l) 

P « M-af/^^+^L (27) 
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Initially, the terms with 5's dominate over the constant terms and the effective equation 
of state of the k-fluid behaves as dark matter: 

a * = ^§^ x0 - <29) 

Later on, as 8 — V 0, the solution tends to the attractor in which the constant terms act as 
a cosmological constant, p A = — M + a^ L ^ 2a l ' . 

On the other hand, the condition (23) states that the speed of sound has to be small 
in order to have the growth in the matter inhomogeneities needed for the structure 
formation. For the field around a minimum the condition is written as 



(30) 



This condition is similar to (26) and, except for very particular Lagrangians, the accom- 
plishment of the first one will be enough. In other words the deviation from the minimum 
5/Xq must be small. 

A different set of Lagrangians that fulfill (22) without the necessity of being around 
the minimum, are those with a large derivative dJf/dX, for example if Jzf = AX 11 , and 
r\ 3> 1. The condition (23) for these powerlaw Lagrangians corresponds to r\ + 1 3> 1 
that is satisfied once the first condition is. 

Another important case is the canonical scalar field Jzf = X — V (0), where the condi- 
tion (22) becomes V((j))/X pa 1. It can be accomplished for different initial conditions of 
the field. The potential V(0) = A/2 [Cosh 2 ( y/3Q/2) + 1] has been studied in Ref. [13] 
due to its property of satisfy this condition. However, the condition (23) corresponding 
to have a small speed of sound is never satisfied as c$ = 1, making it a bad UDM model. 
However see [19] for arguments in favour of the validity of the canonical scalar field as 
UDM. 

CONCLUSIONS 

We have presented a phase space analysis for the unified model (3, 4) of dark matter, 
dark energy, and inflation. We have shown that for a large set of the initial conditions 
(0, 0) a viable dynamics occurs in which inflation (ftty = — 1) happens first, followed by 
a period of dark matter domination (aty = 0), to finish with dark energy (Oty = —1). An 
intermediate radiation period is possible in this model once it is added an extra radiation 
component as in the standard model of cosmology. 

Once inflation ends, the model is fully described by the purely kinetic Lagrangian 
j£f = F (X) with F as in Eq. (3). We have shown that this system possesses a late time 
stable solution in which ftty = — 1, that is dark energy. In ref. [1] the range of parameters 
were given to achieve a successful cosmological model, and in the present work the 
dynamical analysis clearly shows why the system is tenable. A problem however may 
arise in the form of possible instabilities due to the crossing of the system trough 
CQ<p = — 1, as observed in Fig. 1, but this crossing will occur in the future as the equation 
of state should have only moved from aty ~ — 1 in the matter dominated epoch to 
0)(j) ~ —0.75 in the present epoch. 

In the last part of our work we have presented the general features that are necessary 
to have a model that behaves as dark matter. If one adds a cosmological constant to this 
model, one ends with a unified dark matter and dark energy model, called generically 
UDM. There are two conditions that these models should fulfill, equations (22, 23) 
playing the role of an effective fluid with small pressure (in comparison to the density) 
and small speed of sound (in comparison to the speed of light). We have analyzed 
some models studied in the literature that fulfill these conditions. In particular, F(X) 
models that possess a minimum, as Sherrer's model (24) or the model studied in the first 
section (3), when they are close enough to the minimum, they behave as dark matter. 
Departures from the minimum cause a change in the transfer function and therefore 
to a different growth history in comparison to the standard model of cosmology [10]. 
Additionally, we have analyzed other models of the literature, such as models with 
large derivative Jz^, for example ££ = AX 71 with large Tj, and models with canonical 



Lagrangian Jzf = X — V(0). The constrictions were given for these models. 
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